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Abstract 

We discuss skewed parton distributions in the coordinate space. Solution of the 
corresponding LO evolution equation is constructed in terms of eigenfunctions of 
the evolution kernel and its relation to the conformal symmetry is explained. 
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Recently, there has been a renewed interest in QCD evolution of skewed parton distribu- 
tions |1J. Skewed parton distributions play a crucial role in description of hard, exclusive 
QCD processes which are actively considered as tools for investigation of new aspects of 
non-perturbative QCD dynamics. However, it is clear that before non-perturbative infor- 
mation can be reliably extracted from experimental data, all perturbative aspects, such as 
QCD evolution, have to be understood. So far, the main effort has been devoted to studies 
of evolution of the skewed parton distributions (SPD) in the momentum representation 

1- 

SPD are defined through matrix elemets of twist-2 string operators. Consider as an 
example a nonsinglet quark operator: 

0(a,P) = q(^z)zVexp S^-ig z,A^tz)dt^ q(^z) , (1) 

Corresponding SPD can be introduced in the following way 

(P'| 0(a, (3) \P) = N(P') z N(P)e- iar * z J l du F(u, £; /i 2 )e^ (Az) + ... (2) 

where dots denote other Dirac structures. N(P) and N(P') denote initial and final nucleon 
spinors, respectively. The average nucleon momentum is denoted by P = (P + P')/2, 
and the momentum transfer is r = P — P' . The asymmetry parameter £ is defined by 
r • z = 2£P • z. 

Scale dependence of the SPD is governed by a generalised evolution equation 

d n,- + 2\ a s f 1 



fi— F(x, 6*0 = 4^ dyV(x, y, £)F(y, £; /i 2 ) (3) 

At the LO, the evolution kernel V(x, ?/,£) has a set of eigenfuctions associated with local 
conformal operators: 

\cf 2 (x/OV(x,y,Odx = JjCfiy/Z) (4) 

As these eigenfuctions do not form a complete set outside the region > 1, they 

can not be used for expansion of the SPD. This can also be understood by the hybrid 
properties of the SPD. Let us split F(x, £,/x 2 ) in two pieces: 

P(x,£,/i 2 ) = P<(x,£,/i 2 ) + P>(x,£,/i 2 ) (5) 

with 

P < (a:,£, yU 2 ) = 6{x < OP(a;,£,/i 2 ), P>(a;,£,/i 2 ) = 9(x > 0^^^^) (6) 

P<, F > describe partons with x < £ and x > £, respectively. The crucial point is that 
the evolution of P< and F > is qualitatively different Partons which at the initial 

scale belonged to the segment < u < £ stay there in the course of the evolution. On the 
other hand, partons which belonged initially to the segment £ < u < 1 diffuse into the 
segment < u < £ and never come back. Mathematically this means that the function 
F < (x, £, Q 2 ) will be restricted to the initial region < u < £ but P>(x, £, Q 2 ) will expand 
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to whole interwal < u < £. The former and latter cases resemble the ERBL and 
DGLAP evolution, respectively. As it follows, properties of SPD in the region x > £ 
SPD are similar to forward parton distribution f(x), while in the region x < £ SPD looks 
like a distribution amplitude. Expansion in the ortognal set eigenfunctions (Q) is valid 
for the ERBL-region x < £ only, and reflects a typical structure of evolution for such 
configurations of partons. As we see, in momentum space there are two different regions 
in x with different evolution properties. 

The situation is different in coordinate space. Coordinate-space SPD is defined through 
a Fourier transformation: 

^(/?,e;^ 2 ) = -/ 1 ^F(x,e;/i 2 )e^ (7) 

7T J-l 

It is easy to see that, unlike in the momentum space, the coordinate-space SPD J- <t> 
associated with functions (||) are defined in the same interval < f3 < oo. So, one can 
hope that an ortogonal set of coordinate-space eigenfunctions exists and can be used as 
a basis for expansion of SPD. 

Evolution equations in coordinate space have been discussed e.g., in 0, where the role 
of the classical conformal symmetry was emphasized. The authors of || were able to 
write the solution in form of complex integral over conformal spin j. Recently, in [|J the 
solution have been obtained in terms of coordinate-space eigenfuctions corresponding to 
integer j. Here we will obtain solution for the SPD in coordinate space using a formal 
trick and then explain its relation with conformal symmetry. 

Recall that a function f(x) can be expanded in a Neumann series according to || 

00 r°° dX 

f(x) = y £,{2v + 2n)J v+n (x) ^-/(A)J, + „(A). (8) 

n=0 J ° A 

In particular, one finds that e lx/3 can be decomposed according to the Sonine's formula 
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- ) T [3/2] £ i»(3/2 + n) CZ(x/Z) J |+n (^) (9) 

n=0 
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Inserting this expansion in the definition of the coordinate-space skewed quark distribution 
(0) and interchanging summation and integration one obtains: 
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l/2\ 2 °° r 1 s 

f(f3,Z;^)= £' n (3/2 + n)J |+B (/?0 duF(u,^fi 2 )CH^/0- (10) 

V n \KJ „=o Jo 

Now, note that a Gegenbauer moment is proportional to the matrix element of multi- 
plicatively renormalizable local conformal operator and its scale dependence is therefore 
given by 



/* duF(u,£;Q 2 )c!(uj/0 = L l+n /' duF(uj^; fi 2 )d(uj/0, L k 
Jo Jo 



MQ). 



" b o 
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As it follows, the scale dependence of the coordinate space distribution T{J5, £; /i 2 ) is given 
simply by 

3 

.F(/U;Q 2 ) =^ U- ^f(f+n)WJ |+n (^) jf do,F(u;,£; /x 2 ) CS(w/0 .(11) 



n=0 



v 7 ^ V^, 

Now we show that equations (JTTJ) , can be naturally understood as expansions of coordinate- 
space skewed quark distributions in terms of matrix elements of non-local, multiplicatively 
renormalizable, conformal operators. Indeed, applying one can rewrite ([11]) as a 
Neumann-type series: 

F(f3,£;Q 2 )=f3-l ^(3 + 2n)J §+n (/?OW / dXVX ^(X,^,^ 2 ) J, +n (XO (12) 

n=0 " J ° 



We are now in the position to make the relation to conformal symmetry explicit. Let us 
start from the obvious identity 



0(a,/3)= da' df3'6(a-a')5(p - f3')0(a',(3') . (13) 

J-00 Jo 

Applying (H) one finds a representation of a ^-function in terms of a Neumann series 

(35{(3-(3') = J2(l + 2n)A + MJUn(P')- (14) 



n=l 



Inserting this expansion into ([13]) one finds that the string operator 0(a,/3) can be de- 
composed as 

o fOO flic 00 

0(a,/3) = /r* / ^e-^^ + ^J^dmSil^ + j^k;^) (15) 

J— 00 Lll _2 2 

in terms of conformal string operators + j, k; fx 2 ). Such operators 



00 



S(l/2+j,k;^)= j dae tka / d/3 V^WI*^) °(«» ^) > ( 16 ) 



introduced first in [Q, form a representation of a conformal group and are therefore 
multiplicatively renormalizable at a one-loop level 0, i.e. 

5(1/2 + j, fc; Q 2 ) = Lj 5(1/2 + j, fc; /i 2 ) . (17) 

At this point j can be easily identified with the conformal spin. Taking matrix elements 
of both sides of the above equations one immediately reproduces equations ([TT|). 

Note that in the coordinate-space representation the corresponding LO amplitude 
M(£; fi 2 ) can be written in following way: 

POO 

M(£; Q 2 ) ocin d/3 e~ iK T((3^;Q 2 ) (18) 
Jo 

We have checked, using various models of skewed quark distributions, that the numerical 
algorithm for evaluation of physical amplitudes, based on equation (|T8|), gives accurate 
and stable results, see Figure 1 for an example, except for a case where the variable £ 
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becomes small. This is related to the observation that a non-zero £ provides a natural 
cut-off for large (3 behavior of coordinate-space distributions, which significantly improves 
the convergence of the Fourier integral as compared to the forward case. 




u.o u.o ^ 

Figure 1. Typical results of evolution of |M(£)| 2 as a function of £, starting from a £- 
independent initial conditions F(u, £; = l-1641u~2 (1 — -u) 3 - 5 . The solid line denotes 
|M(£)| 2 at the initial scale fi = 1.777 GeV, the dashed line represents |M(£)| 2 evolved to 
fi = 10 GeV. 
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